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The BeppoSAX view of bright Compton-thin Seyfert 2 galaxies 
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Abstract. We present the analysis of 31 observations (17 of which are published here for the first time) of 20 bright 
Compton thin Seyfert 2s, in the 0.1-200 keV band, performed with the BeppoSAX satellite. The sample consists of 
all Seyfert 2s in the BeppoSAX public archive, with a 2-10 keV flux higher than 5xl0 -12 erg cm _2 s _1 . The good 
statistics available and the broad energy band permit a detailed study of the main continuum components of these 
sources, i.e. the primary power-law, the reflected component, the soft emission and the high-energy cut-off. The 
main results of our analysis are: (1) the 3-200 keV intrinsic power-law has a mean photon index F = 1.79 ± 0.01, 
with a dispersion of a = 0.23. (2) The high-energy exponential cut-off at E~ 100 — 300 keV is not an ubiquitous 
property of Seyfert galaxies: in ~ 30% of the objects the continuum power-law does not drop up to energies of 300 
keV or more. (3) A reflected component is present in almost all the sources (17 out of 21). The small variations 
of this component with respect to the intrinsic continuum, in objects with multiple observations, suggests that 
the reflector is not the accretion disk, but must be located much farther from the nucleus. (4) The range of ratios 
between the reflected and intrinsic components suggests that the circumnuclear medium is not homogeneous, and 
a significant fraction of the solid angle is covered by a gas thicker than that along the line of sight. (5) The iron Ka 
line is present is all but one the sources. The equivalent width is in the typical range of Seyfert Is (EW= 100-300 
eV) in sources with low absorption (Nh < 3 x 10 23 cm -2 ), and increases in more absorbed objects, as expected 
according to unified models. (6) The energy resolution of BeppoSAX is in general too low to measure the iron 
line width. However, in 6 cases we measured a significant line broadening. 



Key words. Galaxies: active; Galaxies: Seyfert; X-rays: 
galaxies 

1. Introduction 

In this paper we present the analysis of the observations of 
bright Seyfert 2 galaxies performed with the BeppoSAX 
satellite. The main aim of this work is to investigate 
the average X-ray continuum properties of these sources, 
through a homogeneous analysis of all the BeppoSAX 
archive data. The broad energy band of BeppoSAX (0.1- 
200 keV). permits us to separate the main continuum 
components ( "soft excess" , intrinsic power law, reflection) 
with unprecedented precision. 

Several homogeneous analyses of bright Seyfert 2s 
were published in the past years: Turner & Pounds 1989 
(EXOSAT observations, 0.1-10 keV), Nandra & Pounds 
1994 and Smith & Done 1996 (GINGA observations, 2-20 
keV), Turner et al. 1997 (ASCA observations, 0.5-10 keV). 
Many of the data analyzed in these works have been more 
carefully analyzed using specific models for each source. 
The results are published in many papers where the X- 
ray emission of a single or a couple of sources is studied 
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in detail. However, the results obtained in these detailed 
studies are in most cases in agreement with the simpler 
fits performed in the works quoted above where a large 
number (> 20) of observations are analyzed. The reason 
for this agreement is that even if peculiar features can be 
present in the X-ray spectrum of each source, the statistics 
is dominated by a continuum well fitted by a few compo- 
nents, namely a thermal emission at low energies (from 

0. 1 to a few keV), an absorbed power law with photon 
index r ~ 2, and a reflection component. Therefore, a 
homogeneous analysis of moderately large samples is the 
best way to investigate the average X-ray properties of 
these sources, to test unified models and to look for new 
correlations between the different spectral parameters. 

The main results of the works quoted above support 
the unified models of AGNs, according to which the nar- 
row emission line (type 2) objects differ from broad emis- 
sion line objects (type 1) only for the presence of dusty 
absorbing gas along the line of sight. 

The main X-ray spectral properties of Seyfert 2s can 
be modeled with the following 5 components: 

1. A power-law with photon index r ~ 2. X-ray obser- 
vations of Seyfert Is (Nandra et al. 1997) show that 
this power-law extends down to energies of a frac- 
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tion of kcV. At high energies, an exponential cut-off 
is observed in the BeppoSAX observations of a few 
objects (Guainazzi et al. 1999, Perola et al. 1999, 
2000, Nicastro et al. 2000), at energies from 60 to 
300 keV. The physical origin of this component could 
be a two phase accretion disk, where the soft pho- 
tons emitted by a "cold" (kT < 50 cV) thick disk 
are Comptonized by a hot (kT^ 100 keV) thin corona 
(Haardt & Maraschi 1993). 

2. A photoelectric cut-off of the intrinsic power-law is 
present at some energy, depending of the column den- 
sity, N H , of the absorbing gas. If 10 22 cm~ 2 <N H < 
10 23 cm~ 2 , the cut-off energy is below ~ 2 kcV, if 10 23 
cm~ 2 <Nh < 10 24 cm~ 2 the cut-off energy is between 
2 and 10 keV. If Nh > 10 24 cm _2 s _1 no direct emis- 
sion is observable below 10 keV. The intrinsic com- 
ponent could only be revealed by observations in the 
10-100 kcV band. If N H > 10 25 cm^s^ 1 the source 
is completely thick at any energy, and only reflected 
components (see below) are observable. 

3. Thermal soft component: thermal emission, with kT 
ranging from 0.1 to a few keV is a good description of 
the soft X-ray emission of most Seyfert 2s. This compo- 
nent in Seyfert 2s cannot be due to emission from the 
accretion disk (since otherwise it would be obscured 
like the power-law component), and is probably asso- 
ciated to the warm gas confining the broad emission 
line clouds. 

4. Reflected component: a cold reflected component is of- 
ten required to obtain a good fit of X-ray spectra of 
Seyfert galaxies above ~ 5 keV. The reflecting medium 
could be the accretion disk itself or, alternatively, the 
inner edge of the absorbing gas (Ghisellini et al. 1994), 
or a wind (Elvis 2000). 

A warm reflection component is also required in many 
cases. Assuming that this component is due to scat- 
tering of the primary component with warm electrons, 
the spectral shape is the same as in the incident con- 
tinuum. 

5. Iron line: an iron Ka emission line, at energies E~ 6.4 
keV is observed in most Seyfert galaxies. It is believed 
to be originated both in the accretion disk and in 
the absorbing/reflecting circumnuclear gas (Matt et al. 
1991, 1996) The observed equivalent width of the line 
is typically 100-300 eV in Seyfert Is and in Seyfert 2s 
with Nh lower than a few 10 23 cm~ 2 . In sources with 
heavier absorption the equivalent width is higher (up 
to several keV) since the continuum at the line en- 
ergy is more absorbed than the line, that is -at least 
in part- emitted by the reflecting material along free 
lines of sight. The line component emitted from far re- 
flectors is narrow, and does not immediately follow the 
continuum variations, while the disk component can be 
broad and is expected to vary with almost zero delay 
with respect to the continuum. 

In this work we will use the components described 
above to fit the 0.1-200 keV observations of a sample of 31 



Name 


z 


Class. 


NGC 526a 


0.0295 


Sy 2 


NGC 1365 


0.0055 


Sy 1.8 


IRAS 05189-2524 


0.0426 


Sy 2 


NGC 2110 


0.0076 


Sy 2 


NGC 2992 


0.0077 


Sy 1.9 


MCG-5-23-16 


0.0083 


Sy 2 


NGC 4258 


0.0015 


Sy 1.9 


NGC 4388 


0.0084 


Sy 2 


NGC 4507 


0.0118 


Sy 2 


IRAS 13197-1627 


0.0172 


Sy 1.8 


NGC 5252 


0.0230 


Sy 1.9 


Centaurus A 


0.0018 


Sy 2 


NGC 5506 


0.0062 


Sy 1.9 


NGC 5674 


0.0249 


Sy 1.9 


NGC 6300 


0.0037 


Sy 2 


ESO 103-G35 


0.0133 


Sy 2 


NGC 7172 


0.0087 


Sy 2 


NGC 7314 


0.0047 


Sy 1.9 


NGC 7582 


0.0053 


Sy 2 


NGC 7679 


0.01714 


Sy 2 



Table 1. The sample of Compton thin, X-ray bright 
Seyfert 2 observed with BeppoSAX 



observations of 20 Seyfert 2s, performed by the BeppoSAX 
satellite in the period 1997-2000. 17 out of these 31 ob- 
servations are currently unpublished. In Sect. 2 we review 
the basic properties of the sample and we summarize the 
data reduction process. In Sect. 3 we describe the models 
used to fit the data. In Sect. 4 we discuss the results of 
data analysis. In Sect. 5 we summarize our conclusions. 

2. Data: selection and reduction 

The sample consists of 31 observations of 20 sources (one 
source was observed 5 times, one 3 times and 5 sources 
were observed twice). The basic properties of the sample 
objects are listed in Table 1. 

The only selection criteria are a 2-10 keV measured 
flux higher than 5xl0 -12 erg cm~ 2 s _1 , and a mea- 
sured absorbing column density lower than 10 24 cm~ 2 
(i.e. "Compton thin" sources). As a consequence, the sam- 
ple contains only high signal-to-noise spectra, that can be 
studied in considerable detail. The observation log is re- 
ported in Table 2. All the observations are available in 
the BeppoSAX public archive provided by the ASI Science 
Data Center (SDC). Each source has been observed with 
3 narrow-field instruments on board of the Italian/Dutch 
BeppoSAX satellite: LECS (0.1-10 kcV, Parmar et al. 
1995), MECS (1.65-10.5 kcV, Boella et al. 1995b) and PDS 
(20-200 keV, Frontera et al. 1995). 

The LECS and MECS event files were obtained from 
the SDC public archive. The spectra were created using 
the XSELECT v. 2.0 code. A circular region of radius 4', 
centered on the source, was selected to extract the spec- 
trum. The data were rebinned in order to have at least 20 
counts per bin. This allows the use of Gaussian statistics 
in the model fitting. For the PDS data, we used the spec- 
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Name 


Obs. date 


LECS exp. (s) 


MECS exp. (s) 


PDS exp. (s) 


ReP 


NGC 526a 


31 DEC 98 


38977 


93175 


44199 


1 


NGC 1365 


12 AUG 97 


8862 


27580 


13023 


2 


IRAS 05189-2524 


03 OCT 99 


17378 


41973 


19455 


3 


NGC 2110 


12 OCT 97 


40665 


83701 


38673 


4 


NGC 2992 #1 


01 DEC 97 


26668 


72007 


33820 


5 


NGC 2992 #2 


25 NOV 98 


21914 


59245 


27084 


5 


MCG-5-23-16 


24 APR 98 


35825 


76990 


32524 


UNP. 


NGC 4258 


19 DEC 98 


32954 


99428 


46933 


6 


NGC 4388 #1 


09 JAN 99 


65242 


67878 


51189 


UNP. 


NGC 4388 #2 


03 JAN 00 


20822 


28311 


14607 


UNP. 


NGC 4507 #1 


26 DEC 97 


22229 


55009 


26967 


UNP. 


NGC 4507 #2 


02 JUL 98 


8889 


31408 


16953 


UNP. 


NGC 4507 #3 


13 JAN 99 


11109 


41325 


20081 


UNP. 


IRAS 13197-1627 


22 JUL 98 


14140 


43780 


18616 


UNP. 


Centaurus A #1 


20 FEB 97 


16912 


33393 


14873 


UNP. 


Centaurus A # 2 


06 JAN 98 


16681 


54243 


22987 


UNP. 


Centaurus A # 3 


10 JUL 99 


13676 


38779 


17801 


UNP. 


Centaurus A # 4 


02 AUG 99 


13670 


38756 


17801 


UNP. 


Centaurus A # 5 


08 JAN 00 


11775 


34004 


17941 


UNP. 


NGC 5252 


20 JAN 98 


20897 


61023 


29155 


UNP. 


NGC 5506 #1 


30 JAN 97 


15495 


39363 


17008 


UNP. 


NGC 5506 #2 


14 JAN 98 


8868 


39004 


17500 


UNP. 


NGC 5674 


12 FEB 00 


21029 


45039 


20809 


UNP. 


NGC 6300 


28 AUG 99 


39567 


86265 


38555 


7 


ESO 103-G35 #1 


03 OCT 96 


10088 


50620 


21029 


8 


ESO 103-G35 #2 


14 OCT 97 


3617 


14312 


5915 


8 


NGC 7172 #1 


14 OCT 96 


15212 


39167 


17277 


9 


NGC 7172 #2 


06 NOV 97 


22405 


49185 


21147 


9 


NGC 7314 


08 JUN 99 


35595 


89799 


42564 


UNP. 


NGC 7582 


09 NOV 98 


26859 


56435 


26095 


10 


NGC 7679 


06 DEC 98 


39799 


90908 


45298 


11 



Table 2. a References: UNP: previously unpublished; 1: Landi et al. 2001; 2: Risaliti et al. 2000; 3: Severgnini et al. 
2000; 4: Malaguti et al. 1999; 5: Gilli et al. 2000; 6: Fiore et al. 2001; 7: Guainazzi 2002; 8: Wilkes et al. 2000, 9: Akylas 
et al. 2000; 10: Turner et al. 2000; 11: Delia Ceca et al. 2001. 



trum provided by the SDC public archive. The background 
spectra for the LECS and MECS data were obtained from 
the long observations of blank field provided by the SDC. 
Calibration matrices (redistribution matrix, RMF, for all 
the three instruments, and ancillary (effective area) ma- 
trix, ARF, for LECS and MECS) were also provided by 
the SDC. 

The spectra used in the scientific analysis consist of 
LECS data from 0.1 to 3 keV (above 3 keV the calibration 
of the LECS is uncertain, and we do not lose much in terms 
of statistics, since the area of the MECS is significantly 
higher than that of the LECS above 2 keV); MECS data 
from 1.65 to 10.5 keV, and PDS data from 20 to 200 keV 
(in several cases, where the statistics at high energy was 
too low, we cut the spectrum at 100 keV). 

3. Model fitting 

The data analysis was performed using the XSPEC v. 10.0 
code. All the errors quoted in this works are at a confi- 
dence level of 90% for one interesting parameter. 

In all models we fixed the normalization factor between 
PDS and MECS to 0.8. The uncertainty in this factor is 



at most 10%, according to the SAX SDC. This is not a 
dominant source of error, since te typical final errors in 
the parameters strongly affected by high-energy data are 
typically much higher. This is due to both the limited 
number of counts and the degeneracy between some of the 
spectral components, as will be further discussed in the 
following Subsections. The normalization between LECS 
and MECS were left free. The typical best fit values are 
in the interval 0.6-0.8. 

The Galactic absorption is taken into account in all 
models, by multiplying each component by an extra ab- 
sorption factor. 

First, we fitted the data with a simple model, consist- 
ing of an absorbed power law plus an unabsorbed power 
law. This model is a good zero-order approximation of the 
continuum in Seyfert 2s. Our main aim is to highlight the 
other spectral features from the analysis of the residuals. 
In T able 3 we show the results of the spectral analysis. In 
Fig. k/1 we plot for each spectrum the folded data+model 
and the contribution to x 2 from each spectral bin. This is 
an effective way to look graphically for statistically signif- 
icant unfitted features. 
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Name 


r lux 


1 1 


N H 


1 2 


H 


„ ,2 /j „ r 

X /d.o.t. 


NGC 526a 


1.72 


1 - DD -0.05 


1 c+U.2 
l-6Io.2 


1 o+0. 6 

!-3Io.3 


38 


278/242 


NGC 1365 


0.68 


1 eo+0.21 


44^8 


■r 0+0.3 
i-S-o.s 


11 


129/77 


IRAS 05189-2524 


0.35 


o ^ i +0.26 

2-41Io.3o 


^+2.0 
8-6li.7 


n 1 +0.5 

2 -l-o:e 
1 ^0+0.45 
l-48+ . 25 


218 


88/72 


NGC 2110 


3.00 


1 «e:+0-05 
l.OO_ 04 


n+0.3 
3.9I . 2 


24 


358/282 


NGC 2992 #1 


0.61 


1-46I .08 


n cc+0.13 






253/174 


NGC 2992 #2 


7.37 


-, cr .+0.02 

1-69-o.oa 


0.85I . 04 






402/349 


MCG-5-23-16 


9.28 


J- '^-0.02 


1 ^o+0.05 

1 -63Iq.o5 






536/376 


NGC 4258 


0.79 


9 qo+0.17 

Z.OO_Q 


!3-3+i;i 


1 or+0.10 

1-85+o.n 


23 


-A 1 1 / l rt^7 

266/197 


NGC 4388 #1 


2.53 


1-60Iq:o3 


A -| * + 1.9 


l-73Io:u 


42 


382/239 


NGC 4388 #2 


0.94 


l.t>8_ 012 


541 7 


1 -7C-+0.23 
1-75+0.28 


18 


166/105 


NGC 4507 #1 


1.84 


1 co+0.05 
1-68I .0B 


64I 4 


1 9+0-2 
i - Z -0.2 


174 


334/186 


NGC 4507 #2 


1.65 


1 cc+0.07 

1-66I . 08 


rn+5 

59_5 


n OQ+0.35 

u.y^_ 25 


260 


230/138 


NGC 4507 #3 


0.87 


l-SOIo.io 


641 10 


1 1+0.3 

1-1-0.3 


63 


258/127 


IRAS 13197-1627 


0.47 


1 vn+0.33 

1-70I . 28 


OQ + 7 


q c+1-2 
J '°-0.9 


6.5 


1 A A 1 Cl T 

144/87 


Centaurus A #1 


19.6 


1 7Q+001 
1 - 1 y -0.02 


1 n 1 +0.3 
10-1-0.2 


1 ec+0.23 
l-8O_ 24 


56 


417/259 


Centaurus A #2 


25.4 


1 r,r, + 0.01 


n /in+0.16 

9-40I 015 


9 nc - + 0.23 
Z.UO_ 25 


67 


396/273 


Centaurus A #3 


23.6 


1 „ o +0.01 

1-78Iq.oi 


n qo+0.21 

y-o^-o^o 


1.89^0.2? 


65 


403/255 


Centaurus A #4 


23.6 


1 „ o +0.01 

l-78I .oi 


n oc+0.19 

9.86I 019 


9 i n +0.20 
^• 1U -0.23 


54 


407/266 


Centaurus A #5 


23.1 


1 „„+o.02 


n no+0.24 

9.93I .23 


1 cq+0.23 


70 


374/247 


NGC 5252 


0.26 


-i ro+030 
i-O^-O.SO 


e o+1.8 
6-2IL3 


3-7_2.4 


29 


80/70 


NGC 5506 #1 


7.47 


1 QO +0.05 
t.8^_o.04 


3-46Io. 15 


1 n+0.3 

l-2Io.2 


96 


408/256 


NGC 5506 #2 


7.21 


1 7K+0-03 
^^-O.OS 


9 11 +0.16 

3-41I .i8 


1 -7+4.0 

1-7-o.s 


76 


340/228 


tvt/^i/~i pit/"* ^ 

NGC 5764 


0.54 


1 - n +0.21 

1 -70Io.i9 


^ + 1.7 

6-8Ii.3 


o n+0.9 

2-9Ii.o 


13 


91/107 






-1 7-1+0.07 
± -' ± -0.07 


01 O + l -2 

zi.y x 2 


9 0+0.3 

2 -°-0.3 


7°. 


ZDl/ loo 


ESO 103-G35 #1 


2.63 


1 O4+O.O6 


20.5+°j 






255/183 


ESO 103-G35 #2 


1.53 


-1 70+0. 18 
± -'°-0.14 


20+* 






56/76 


NGC 7172 #1 


1.12 


1.68l° ;li 


10.5+1 


2 6+°' 7 

z -°-0.6 


123 


152/165 


NGC 7172 #2 


0.60 




s.9±i:S 


9 q+0.5 
z -°-0.7 


46 


153/161 


NGC 7314 


2.39 


1 04+0.07 


i-08l° ;^ 


1 n3+ 007 

1 - uo -0.09 


25 


339/288 


NGC 7582 


2.0 


1 cc+0.05 
1 -°°-0.05 


12.8«; 8 8 


9 1+0.3 
z - x -0.3 


36 


265/187 


NGC 7679 


0.59 


1 q^+0. 23 

i.ao_ 20 


o.io±S°I 


1 0+°' 3 


14 


207/222 



Table 3. Model A: Absorbed power-law plus a second unabsorbed power-law. Since this second component fits mainly the 
soft excess, the ratio between the normalizations of the two components gives an idea of the relative strength of the soft excess 
with respect to the primary component. a Flux 2-10 keV in units of 10 11 erg s _1 cm -2 ; b Fi: Photon index of the absorbed 
power-law; c Nh: X-ray absorbing column density, in units of 10 22 cm -2 ; d F2'- Photon index of the transmitted power-law; e 
R: Ratio between the normalizations of the absorbed and transmitted components. 



From the analysis of Fig. |A.2| we can identify the fol- 
lowing unfitted features: 

— At low energies (E<2 keV) several sources present a 
more complex spectral shape than a simple power- 
law. This is apparent, for example, for the objects 
NGC 526a, NGC 4258, NGC 4388, NGC 4507, Cen A. 

— An excess around 6.5 keV is present in most observa- 
tions. 

— The high energy spectrum (E> 20 keV) is in many 
cases convex, and cannot be reproduced by a single 
power-law. 

We also note that in all the sources but one (NGC 7679) 
the photoelectric cut-off is highly significant, and permits 
a precise determination of the absorbing column density, 
N H . 

At high energies the observed curvature can have two 
different origins: an high energy exponential cut-off, or a 
reflection component. In this case a degeneracy between 
these two components is possible, so we fitted our data 
with two new models: in the first one (model B) we re- 



placed the unabsorbed power-law with a Raymond-Smith 
thermal model, and we added a Gaussian and an expo- 
nential cut-off to the absorbed power-law. In the second 
model (model C) we also added a warm reflection compo- 
nent, using a power-law with the photon index equal to 
that of the primary continuum, and a cold reflection com- 
ponent, using the PEXRAV model in XSPEC. This model 
assumes an infinite plane geometry for the reflector. This 
may not be a good approximation of the real reflector. 
However, the overall spectral shape of the reflection com- 
ponent should be only weakly dependent on the details of 
the absorber. The dependence of the spectral shape on the 
reflection angle in the PEXRAV model is low, therefore, 
given the uncertainty in the real absorber geometry, we 
chose to fix the inclination angle in the PEXRAV model 
to the value of 9 =30 deg (the value has been chosen in 
order to have an average ratio between the observed sec- 
tion of the absorber and it real area of cos6> = 0.5). The 
results of the fits with model C are summarized in Table 
5. 
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Name T a N H E c c kT d E| e Wj c EW| C X 2 /d.o.f. 



NGC 526a 


1 cu +u,ub 
i - o *-0.04 




> 176 


>1.1 


6.65!-- 


n 00 + 21 
U.OO_ 33 


1821^ 


251/238 


NGC 1365 


() 7 r+0.23 
u -'°-0.45 




44lf 


cy 7 + 3.O 
Z ' ' -1.1 


6.2212:12 


- 


642+ 248 

0*^-230 


106/74 


IRAS 05189-2524 


9 OO + 0.26 
Z.OO_ Q 30 


7 6 +15 


> 15 


0.881°;^ 


fi 47+0.23 
D - 4 ' -0.20 


< 0.45 


oon+187 
^ oa _153 


86/68 


NGC 2110 


1 62+ 07 


7+0.3 
°- ' -0.2 


> 166 


> 2.6 


n q^l+0.07 

o-J^-o.os 


< 0.32 


9qc+63 
Z JJ_60 


275/278 


NGC 2992 #1 


1 RS+ 008 


i-i±8:5 


> 112 


1 1+ 1 - 2 

J - J '-Q.2 


6.611^7 


1fi +0 11 


684111 


172/170 


NGC 2992 #2 


1 71 +0-02 
' ^-0.04 


0-851°;^ 


> 215 


- 


6.551^1 


m 


10311* 


359/346 


MCG-5-23-16 


1 RQ+ 003 


1 C4+0.07 


26lli 2 


- 


6.44lg° 6 6 


— 


108111 


440/373 


NGC 4258 


9 9c:+0.03 
^-^-o.os 


i2.i±g:5 


> 145 


67 +0 - 74 
u -°' -0.07 


6-56l2:? 9 


< 0.35 


80tg 


162/192 


NGC 4388 #1 


1 90 + O.O8 
± - zo -0.10 


31.4±1;° 


108+ 42 


9 a+0.6 
Z '*-0.6 


a 49+O.O7 
"•^-0.07 


0.2812:^ 


606±1| 4 


272/236 


NGC 4388 #2 


1 9^+0.18 
1 - ZJ -0.11 


38l 4 3 


264 +«° 


4+2.2 


6 S.6+ 0,05 
°- JO -0.05 


< 0.12 


5641H2 
11501H2 


111/101 


NGC 4507 #1 


o-69±g:i| 


38.5 +4 : 8 5 


36l 9 8 


3.7±f$ 


O-JJ-0.10 


o.43i2:ll 


234/182 


NGC 4507 #2 


0.56±°;g 


34.61J-2 


35l 3 3 


> 2.4 


fi 44+ 08 
D - 44 -0.09 


0.4412.1? 


940lH 8 3 


150/134 


NGC 4507 #3 




20.41™ 


27.i±i:i 


Q+4.5 
z -°-1.0 


fi 44+ ' 06 
D -**-0.06 


fl 90+0 09 
u - ZJ -0.08 


loooll™ 


164/124 


IRAS 13197-1627 


1 qi +0.08 
L - O1 -0.23 


24.4l^j 


> 106 


LOIS! 


fi 4n+ 09 

D - w -0.08 


m 


463llf 


116/77 


Centaurus A #1 


1 79+0.03 
< ^-0.04 


y - ' u -0.26 


> 231 


6-71^ 


6-4ol2:2 6 6 




98+1? 


297/256 


Centaurus A #2 


1 7Q+0.02 
1 -' u -0.03 


q Q7+0.17 

y ' ul -0.09 


351+^1° 


4-7l?;l 




o m 


50+11 


342/270 


Centaurus A #3 


i 70+O.O2 
1 -'°-0.03 


9 -^°-0.22 


> 336 


5.9±l:§ 


c C9+O.O4 
°- oz -0.05 


o m 


109111 


271/252 


Centaurus A #4 


1 7O+0.03 
1 -'°-0.03 


q K4+0-22 

9 - J ^-0.23 


> 345 


4.oli;l 


c C9+O.O5 
"•0^-0.04 




10811* 


291/263 


Centaurus A #5 


1 QO+0.02 


10 1+ ' 2 


> 588 


15l^ 2 


fi 4R+ 05 
D - 4D -0.05 


m 


H8112 


307/244 


NGC 5252 




6 1+ 1 - 4 


> 58 


> 8.6 


6.48±g:H 


< 0.24 


267±1|? 


76/66 


NGC 5506 #1 


1 09+0.02 
1 ' oz -0.03 


90+O.13 
°- io -0.12 


> 415 


o.5oi°: 6 2 ; 


6.59±g:°| 


9R+ 012 
u - zo -0.09 


199l 5 3 2 


321/253 


NGC 5506 #2 




O 1 r + U.14 

°- lo -0.21 


214l^f 


9 q+4.4 
— 2.2 


6 48+ 0,08 


< 0.22 


130l 4 ^ 


282/224 


NGC 5764 


l-69Io.i4 


6-3to. 7 


> 340 


-* -i +2 3 
1 1 T^-o 

J— 1 -0.3 


6.36l2ig 


< 0.47 


zz± -108 


82/103 


NGC 6300 


1 62+ 07 


20.7+1? 


> 184 


1 0+ 17 


6.4 m 


44+ - 26 


267ll|? 


236/180 


ESO 103-G35 #1 


1 69+° 13 


18.2+? : « 


84lJ° 4 




6.44l2:f 2 


20+ - 20 

u - zu -0.10 


234+ 424 


185/181 


ESO 103-G35 #2 


1 62+° 19 

± - uz -0.32 


18.3+11 


> 67 


> 0.22 


6.49l2:l 2 9 


< 0.35 


244lll 9 5 


44/72 


NGC 7172 #1 


1 64+° 12 


q 7 +o.9 

y -'-0.8 


> 109 


< 1.1 


6.4612:" 


< 0.6 


1161 427 


148/162 


NGC 7172 #2 


1 79+O.O8 
'-■'^-O.ll 


10.8+?'? 


> 103 


> 2.4 


6 62+ 09 


< 0.33 


1571 4 ? 1 


149/157 


NGC 7314 


1 7 n+o.05 

3 -0.03 


0.86«;^ 


> 504 


n 10+0. 05 
0-18_ .03 


6.49l2:l 3 3 


< 0.42 


169+^ 
1121 4 ? 


330/284 


NGC 7582 


-1 01 +003 
± - J± -0.03 


io.5t$:* 


10911° 


1 04+ ' 08 


6.2812:1? 


m 


243/184 


NGC 7679 


1 7n+ 007 

1 - lu -0.10 


< 0.05 


> 90 


1 1+0.9 
ul -0.2 








208/221 



Table 4. Model B: Absorbed power-law, Raymond-Smith component, iron line. a F: Photon index of the absorbed power-law; 
b Nh: X-ray absorbing column density, in units of 10 22 cm -2 ; c Ec: Exponential cut-off energy (keV); d kT: Temperature of the 
thermal component in units of kT (keV); e Ef c : Fe Ka line energy (keV); ^ Wf c : Line width (keV); 3 EWfb: Line equivalent 
width (eV). m Fixed parameter. 



4. Results and Discussion 

Our models A,B,C provide a good fit to almost all the 
sources of the sample: we obtained a reduced % 2 lower 
than 1.10 in 29 cases out of 31. The best model for each 
observation, and the reduced x 2 are listed in Table 6. We 
briefly comment about the "bad" fits (x 2 > 1-1) in the 
Appendix. Here we discuss the main results of our analy- 
sis: 

4.1. Photon Index 

The wide energy range and the good statistics of our spec- 
tra allow a precise determination of the photon index of 
the continuum power-law: 90% errors are typically a few 
percent. The photon index distribution is plotted in Fig. 
[l] versus the ratio between the reflected and transmitted 
component (see below for a discussion on this ratio). The 
data are from the best fit model, either B or C, according 
to the minimun reduced % 2 . The average photon index is 
T A v = 1-79 ± 0.01. The la dispersion is 0.23. A similarly 
interesting parameter is the "effective" photon index, i.e. 



the slope of a simple power-law fitting the whole spec- 
trum (model A). This parameter is important for studies 
where the average X-ray SED of AGNs is approximated 
with a simple power-law (for example, in most synthesis 
models of the cosmic background). The average photon 
index derived from Table 3 is Te = 1.76 ± 0.01 with a 
la dispersion of 0.21. The small discrepancy between the 
two evaluations of T are due mainly to the contribution 
of the reflection component, which is higher at high ener- 
gies, so flattening the spectrum. The contribution of the 
high-energy cut-off is small, since for most of the sources 
only a lower limit can be determined (see below). We note 
that our sample is the best currently available to deter- 
mine the average T of Seyfert galaxies. The typical extent 
of this interval is from 2-3 keV up to 200 keV, and even in 
the worst cases it is more than a decade (from 5-6 keV to 
100 keV for NGC 4507). This is a significant improvement 
with respect to past works, based on GINGA or ASCA 
data. The photoelectric cut-off does not significantly re- 
duce the energy range over which we observe the intrinsic 
emission of the AGN and, on the other side, makes simpler 
the separation of the low energy component. 
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Name 




1 


N H 


liC 


k 1 






\H if 

W Fo 


EW Fo 


r> h 

Kc 


r> h 


X /d.o.l. 


NGC 526a 


1 


r 4+11.00 

•54I .o4 


L6-0.2 


> 176 


>1.1 


6. 


cc+O.lV 
OO_ 21 


n 00+0.21 

U.>M_ 33 


182_ 96 


< 0.5 


< 0.1 


OKI / O O >T 

251/237 


NGC 1365 


1 


-7,1+0.21 
• 74 -0.39 


oq+20 


> 56 


1 r\+0.6 
1-0-0.4 


6 


n^+0.11 

■24+ . 12 





onn+129 

290I 139 


3-8l 2 


< 0.1 


91/72 


TT» AO O 1 — I OA firn ,\ 

IRAb 05189-2524 


2 


cq + 0.35 
Oy -0.43 


8 - 6 -o.7 


> 25 




6 


'46I .25 


< 0.45 


occ+210 

256 Ii92 


< 4 


n r\o+0.02 

O-02+o.oi 


78/68 


JNGC 2110 


1 


■ 66 -o.o9 


Q Q + 0.3 

d -°-0.3 


> 143 


> 2.9 


6 


oo+0. 08 

■^"-o.os 


< 0.33 


O A + 6 1 


< 0.4 


< 0.03 


274/277 


TvT/~i/~i oooo ii-i 

NGC 2992 #1 


1 


77+O.I8 
' ' ' -0.16 


1 a+0A 
!- 4 -0.4 


> 53 


1 1/^ + 1.06 
l-16-Io.29 


6 


^q+0-06 
D,3 -0.08 


< 0.28 


r^^r + 160 

765+ 8 g u 


n o+4 

O-8+.o.e 


< 0.06 


170/169 


"\T/~i /~i o o o o / / o 

JNGC 2992 #2 


1 


7O+0.08 
• ' y -0.03 




> 141 




6 


c-9+0.14 
OZ -0.33 


n oc+0.20 

0.36+Q12 


1 c/i+47 
154I 33 


n o+0.4 

0-3I .i 


< 0.01 


O f ft 1 O A A 

350/344 


MCG-5-23-16 


1 


77+0.04 
' ' ' -0.04 


1 e;Q,+ - 04 

i.oy_ g 4 


1571^4 




6 


a r+0.08 

'45I .o7 


< 0.30 


1 0/^+31 
126+^ 


r\ r ^+0 2(i 

0-54+^.i 7 


< 0.01 


A f\ A 1 O ^70 

404/373 


JNGC 4258 


2 


OC+0.04 
' /O -0.05 


1 1 +0.7 
12-1+o.e 


> 145 


n c-7+0.07 

0-67-0.07 


6 


r R +0.30 
OD -0.28 


< 0.35 


on+70 

8O+40 


< 0.6 


< 0.02 


162/191 


JNGC 4388 #1 


1 


ro+0.08 
' OO -0.22 


38+4 


> 145 


1 ^n+0.25 


6 


•46I .io 


< 0.23 


900 + 115 


A O+0 2 

0-3l .i 


n m c+0.004 
0.016I .003 


236/233 


"\T/~i /~i A O O O / / O 

JNGC 4388 #2 


1 


a 7+0.04 

• 4 7I .4i 


A 0+I8 


> 109 


1 (T7+0.28 
1-07I .22 


6 


00+0.05 
■^"-o.oe 


< 0.12 


r r + 115 

525tn2 


< 1.4 


n nc+0.02 

0.06+0.02 


100/99 


JNGC 4507 #1 


1 


m +0.06 
' O1 -0.05 


cn+8 


1 oe+153 
126+4 8 


A C7+0.20 

0.87+ 014 


6 


^n+0.10 
•40+0.19 


n oo+0. 14 

0-33+o.ie 


on -7+109 
• 1(1/ 1 

z,£ " -100 


n -7+0.2 

0-7+0.2 


< 0.01 


OA/ 1 / 1 OO 

206/180 


"\T/~i /~1 /I f OT II o 

JNGC 4507 #2 


1 


rc+0.05 


541^ 


> 100 


n -71 +0.27 

0-71+o.4o 


6 


^-7 + 0.11 

•47+0.12 


n 07+0. 18 

0-37+o.i 6 


520lf^ 


n n+0.3 

0-9I .3 


< 0.01 


1 OO / 1 OO 

132/133 


NGC 4507 #3 


1 


-7n+0.04 

•70I .i3 


T1 +20 
71-16 


> 216 


U - '"-0.26 


6 




< 0.25 


a oo + 242 

433I 29 


n r>+0.5 

2-0+o.s 


< 0.02 


1 00 / 1 00 

132/122 


IRAS 13197-1627 


1 


77+0.07 
' ' ' -0.07 


33^4 


> 164 


n nn+0.13 

0-90+o.io 


6 


^0+0.10 


< 0.27 


roc+275 
oeo_i9o 


9 1 +0.5 


< 0.01 


98/77 


Centaurus A #1 


1 


01 +0.05 
• 8i -0.03 


1 n n+0.3 
10-0-0.2 


> 429 




6 


^o+0.08 

•48+0.07 


n oo+0.12 

0-33Io.ii 


1 -70+39 
173+30 


n 1 O+0.19 

0-18+0.07 


n ni r+0 002 
u.uio_ 002 


f / f j 

259/254 


Centaurus A #2 


1 


7O+0.03 
■' 8 -0.02 


Q on+0.13 


> 490 


n cc+0.15 

0.56I . 2g 


6 


+0.09 
•41-0.09 


< 0.5 


-7K+36 
75+27 


n 1 c+0.07 

0-15+0.05 


n ni i +0.002 
0-011I .002 


OOO / O/" 1 T 

290/267 


Centaurus A #3 


1 


70+O.OI 
• ' d -0.01 


O ^0 — \~ ().2i5 

9-43-o.os 


> 339 




6 


t-7+0.05 
-0.04 


< 0.14 


1 00+14 

138+ 21 


< 0.07 


n ni c+0.002 

U.U10_o 002 


f / f 

259/250 


Centaurus A #4 


1 


70+0.02 
■ ■ ^-0.01 




rqA+205 

0'3U_ lg6 


n c/i+0.15 
0- 54 I .26 


6 


C-7+O.O2 
°' -0.08 


< 0.19 


1 /ti +17 

141-21 


< 0.08 


n ni zi+0.002 

0-014I 002 


245/260 


Centaurus A #5 


1 


7C+0.01 
• '°-0.01 


n /IO+0.14 

9.43+ on 


> 605 


0-47I .24 


6 


A r+0.04 

•45+o.oe 


< 0.11 


1 1 T+20 
H7-16 


n nn+0.02 

0-09+0.02 


n ni a +0.002 

0-014+0.Q02 


orr / /I 1 

255/241 


a t r~\ f or o 

NGC 5252 


1 


go+0.13 
• 8d -0.12 


n + 1.6 

6.8+ ? 


> 49 


n e a +0.47 
0-54^0.34 


6 


A T+0. 15 

•47Io.i7 


< 0.25 




1 n + 1.2 

I.6+1.2 


< 0.01 


71/65 


NGC 5506 #1 


2 


no+0.02 
,U,3 -0.03 


/?o+0.17 

■J-o^-o.og 


> 398 


> 0.27 


6 


7/^+0.09 
' U -0.14 


r> O/- + 0.15 

0.36+^i° 


n n r\ — \- H 2 

2 29I 58 


O O + l 


< 0.01 


269/251 


NGC 5506 #2 


2 


no+0.02 
,U ^-0.02 


oc;+0.10 
,; '- 80 -0.19 


> 314 


n ,1+1.0 

0-84+o.s 


6 


C9+0.15 
OZ -0.11 


< 0.25 


1 n j+103 

164_ 3 2 


-t 0+0.2 
l-2Io.2 


< 0.03 


O J /ooo 

224/220 


JNGC 5674 


1 


oo+0. 09 
■° -0.19 


a c+0.8 

6-5I .5 


> 112 


1 1+1-7 

l-!-o.7 


6 


oe+0.17 
•36I .37 


< 0.5 


oor+219 

236+ 124 


1 + 2 

-0.7 


< 0.06 


TO / 1 OO 

79/102 


IN^r^ DoUU 


i 


07+O.O4 
•°' -0.04 


99 0+0.9 

-"•°_0 9 




n 09+0.07 

0.8/_ .09 


D. 


97+0.12 
z 1 -0.12 


u 


1 c;n+ 45 

10U_44 


1 1+ 02 
1 - i -0.2 


>• n m 

U.U1 


1 Q7 / I 70, 

iy 1 / ± / y 


ESO 103-G35 #1 


1 


oe+0.09 
•°°-0.18 


19.7t° : ? 


87l 33 4 




6 


07+O.13 
•°' -0.09 


22+ ' 16 

u - zz -0.05 


215l 6 9 2 


4+ ' 2 


< 0.01 


172/180 


ESO 103-G35 #2 


1 


oc+0.13 
° J -0.33 


20.5j£J 


> 84 




6 


Ziq+0.12 
'^ 3 -0.15 


< 0.33 


224U« 


7+ ' 8 
u - ' -0.6 


< 0.01 


42/73 


NGC 7172 #1 


1 


oo+0. 06 
' 88 -0.10 


11 1+ 12 

2 


> 128 


< 0.44 




6.4 m 


m 


901™ 


1 1+ 2 


< 0.01 


142/162 


NGC 7172 #2 


1 


Q7+0.09 
a ' -0.08 


11 0+ ' 9 


> 105 


1Q+ 004 

u - ±a -0.08 


6 


rr+0.17 
JJ -0.13 


< 0.31 


199litt 


9 4+0.8 

z -^-0.9 


< 0.01 


140/156 


NGC 7314 


2 


nv+0.02 

■ U ' -0.05 


1 99+O.O9 
1 - zz -0.14 


> 462 


7+4.0 
°- ' -1.3 


6 


07+0. 18 
°' -0.17 


< 0.45 


102 + |f 


9 4 +o.4 
z -^-i.o 


< 0.04 


276/283 


NGC 7582 


1 


QO+0.03 
' ao -0.03 


14 3+ a7 


> 364 


0.78+°° 9 2 




6.4 m 


m 


1831^ 


9 9+O.6 
z - z -0.2 


< 0.01 


175/184 


NGC 7679 


2 


qo+0.24 
u °-0.21 


6+ 05 

U.O_ .4 


> 38 












2-7±?:| 




198/220 



Table 5. Model C: Absorbed power-law, Raymond Smith component, cold reflection, warm reflection, iron line. a V: Photon 
index of the absorbed power-law; 6 Nh: X-ray absorbing column density, in units of 10 21 cm~ 2 ; c Eq: Exponential cut-off energy 
(keV); d kT: Temperature of the thermal component in units of kT (keV); e Ef c : Fe Ka line energy (keV); f Wf c : Line width 
(keV); 9 EWf g : Line equivalent width (eV); h Rc, Rw: Ratio between the normalizations of the reflected and direct components. 
m Fixed parameter. 



4.2. Reflection components 

Cold reflection: the cold reflection component peaks 
around 30 keV. Therefore, the PDS instrument allows a 
much better measurement of this feature than low energy 
instruments. The model we used assumes a viewing angle 
of 30°. Given the statistics in the PDS range, a variation 
of this angle is analogous to a change in the normaliza- 
tion, since the spectral changes are too tiny to be dis- 
criminated. The quantity used to measure the strength of 
reflection is the ratio, R, between the normalizations of 
the reflected and transmitted components. A comparison 
between Table 4 (model B, with no reflection) and Table 5 
(model C, with reflection) illustrates the importance of a 
correct estimate of the reflection component. In most cases 
adding a reflection component significantly improves the 
fit (26 out of 31, Table 6). Moreover, in these cases a fit 
without reflection gives biased estimates of the continuum 
components: for example, photon indexes and high energy 
cut-offs are systematically lower in Table 4 than in Table 
5. 



Finally, it is worth noting that the errors associated 
to R are significantly higher than those of the other free 
parameters. This is an effect of the partial degeneracy be- 
tween the different spectral components, in particular the 
intrinsic power law, the high-energy cutoff and the cold 
reflection component. If we freeze these parameters when 
estimating the errors on R, we obtain much smaller errors. 

One of the main debates regarding the cold reflection 
component is whether it originates from the accretion disk 
or from cold material farther from the center (eg. the puta- 
tive "torus" of unified models). Our sample contains mul- 
tiple observations of 7 sources. This gives us the possibil- 
ity to study the long term variability of the cold reflection 
component, and thus to give an estimate to the distance of 
the reflecting medium: if the reflector is the accretion disk, 
we expect the reflection and transmitted component to be 
closely related. In the other case, if the reflector distance 
from the center is greater than the light path in the char- 
acteristic variability time of the primary component, we 
expect the reflected component to remain constant. In or- 
der to study the correlation between the variability of the 



G. Risaliti: BeppoSAX observations of Seyfert 2s 



7 



Name 



Mod. 



NGC 526a 


B 


1 


05 


NGC 1365 


C 


1 


26 


IRAS 05189-2524 


C 


1 


15 


NGC 2110 


B 





99 


NGC 2992 #1 


B 


1 


01 


NGC 2992 #2 


C 


1 


02 


MCG-5-23-16 


C 


1 


08 


NGC 4258 


B 





84 


NGC 4388 #1 


C 


1 


01 


NGC 4388 #2 


C 


1 


01 


NGC 4507 #1 


c 


1 


08 


NGC 4507 #2 


c 





99 


NGC 4507 #3 


c 


1 


08 


IRAS 13197-1627 


c 


1 


27 


Centaurus A #1 


c 


1 


02 


Centaurus A # 2 


c 


1 


09 


Centaurus A # 3 


c 


1 


04 


Centaurus A # 4 


c 





94 


Centaurus A # 5 


c 


1 


10 


NGC 5252 


c 


1 


09 


NGC 5506 #1 


c 


1 


07 


NGC 5506 #2 


c 


1 


02 


NGC 5674 


c 





77 


NGC 6300 


c 


1 


10 


ESO 103-G35 #1 


c 





95 


ESO 103-G35 #2 


B 





61 


NGC 7172 #1 


c 





88 


NGC 7172 #2 


c 





90 


NGC 7314 


c 





98 


NGC 7582 


c 





95 


NGC 7679 


c 





90 



Table 6. Summary of best fit results. Column 2: best fit 
model; column 3: reduced \ 2 f° r the best fit model. 



transmitted and reflected components, we define for each 
pair of observations of the same source a ratio f=Fi/F2 
where F is the 2-10 keV flux, and a ratio r=R.2/Ri. If the 
reflected component varies together with the transmitted 
one, the quantities Ri are expected to remain constant, 
and therefore r=l. If, on the other side, the absolute flux 
of the reflected component remains constant, then R; oc 
1/Fj, and r=f. In Fig. [2] we plot the quantities r and f. For 
each pair, we chose the fluxes in order to have F2 < F\ 
(then f> 1). For the two sources with more than two ob- 
servations (NGC 4507, three observations, and Cen A, 5 
observations), we have 2 and 4 pairs, respectively. In these 
cases, F2 is defined as the lowest measured flux among all 
the observations of the same source. We also plot the two 
lines corresponding to r=l and r=f. The plot clearly shows 
that the r=f line better represents the measures. This is 
strong evidence supporting a distant reflector, while the 
cold reflection from the disk is ruled out by the data. 

If the cold reflection is mainly due to a distant reflec- 
tor, the simplest structure of the circumnuclear medium 
would be obtained with a single cold medium responsible 
for both the cold reflection and absorption of the hard 
X-rays. However, in this scenario the ratios between the 
reflected and transmitted components are somewhat sur- 



2.5 



2 - 




Fig. 1. Photon index versus the ratio R between the normal- 
izations of the reflected and transmitted components. The line 
represents the average value of T (Fav = 1-80 ± 0.03). 



prising. The model we used assumes a perfect efficiency 
reflection (i.e. no transmission) from an infinite plane slab 
covering a 2tt angle. The normalization factor of the model 
is relative to the intrinsic continuum. Therefore, even as- 
suming a perfectly face-on line of sight, the ratio be- 
tween the normalizations of the reflection model and of 
the power-law should be not higher than 2. The value 
R=2 is for a completely Compton-thick reflector cover- 
ing the whole plane angle. Assuming a different geometry 
for the reflector, in particular a non-planar one, can only 
slightly change the overall efficiency, as we can evaluate 
from the comparison of the PEXRAV model with the one 
of Ghisellini et al. (1994), which assumes a toroidal geom- 
etry for the absorber: the results are in agreement within 
a factor of ~ 2. The model of Ghisellini et al. 1994 also 
shows that the reflection efficiency drops with the thick- 
ness of the reflector: at 30 keV the reflection efficiency is 
55% for a Compton thick reflector (N H ~ 10 25 cm -2 ), 
23% for Nh ~ 10 24 cm -2 and only 8% for N H ~ 10 23 
cm -2 . The values for the ratio we obtained (Table 5) are 
in many cases too high to be explained assuming a ho- 
mogeneous reflector having the same column density as 
measured in absorption. In the extreme case of NGC 1365 
the best fit value, R=3.8, is not acceptable. We discuss 
this case in the Appendix. Out of the other 20 objects, 
17 have the value of R higher than the maximum allowed 
by the models described above, the only exceptions being 
NGC 526a, NGC 2110, NGC 4258 and Cen A. Part of the 
discrepancy can be due, in some cases, to the fact that the 
source is observed in a low state, and therefore the value 
of R is high, as explained above. However, this cannot be 
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F./F 

Fig. 2. Reflected component variations versus flux variations 
(see text for definitions) The solid line is the correlation ex- 
pected if the cold reflection component remains constant when 
the intrinsic continuum varies. The dashed line is the corre- 
lation expected if a variation in the intrinsic emission is im- 
mediately followed by a subsequent variation of the reflected 
component. Sources with more than two observations are plot- 
ted as triangles (Cen A) and squares (NGC 4507). 



the case for 18 sources out of 21. Moreover, even in those 
cases where multiple observations are available, R is too 
high even in the high flux states. The only possible ex- 
planation is that the absorber is not homogeneous, and a 
Compton thick medium covers a significant fraction of the 
solid angle (but not the line of sight). We note that other 
recent studies on variability of Nh support the view of a 
non homogeneous absorber (Risaliti et al. 2002). 

Finally, our sample is useful to investigate the puta- 
tive correlation between the photon index and the ratio R. 
Such correlation has been first claimed by Zdziarski et al. 

1999, and further analyzed by Matt (2000) and Petrucci et 
al. (2001), taking advantage of the broad BeppoSAX band. 
At present it is still debated, since the two parameters are 
not independent in the fitting procedure, an increase of 
r being compensated by a decrease of R. In Fig. [j] we 
plot the two quantities. A slight increase of T for R> 2 is 
present, but the correlation is weaker than that found by 
Petrucci et al. 2001 for a sample of bright Seyfert Is. We 
can explain this discrepancy using an argument of Matt 

2000. In order to give a warning about the reality of this 
correlation, this author shows that if the soft component 
of the Seyfert 1 MKN 841 is fitted with a warm absorber 
instead of a blackbody, the best fit photon index changes 
by ~ 0.4 and, as a consequence, R changes by a factor or 
~ 2. This example shows how strong is the dependence 
among the different components of the continuum mod- 
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Fig. 3. High energy cut-off versus photon index of the primary 
power law. 



els. As we already pointed out above, our fits are made 
safer by the presence of the photoelectric cut-off, which 
clearly separates the soft component from the others. As 
a consequence, the statistical dependence among the two 
quantities is weaker in our sample. Therefore, the results 
plotted in Fig. |l| suggest that the correlation is probably 
not real. 

Warm reflection: Adding a warm reflection compo- 
nent improves the fit only for three sources (IRAS 05189- 
2524, NGC 4388 and Cen A, see Table 5). In these cases 
the reflection efficiency is between 1% and 6%. In all the 
other cases we were able to estimate an upper limit for 
the warm reflection of the order of a few percent of the 
intrinsic component. 



4.3. The high-energy cut-off 

The energy of the exponential cut-off is plotted versus the 
photon index in Fig. ||. In most cases we have only lower 
limits, however some interesting conclusions can never- 
theless be drawn. Many lower limits are too low to give 
significant constraints. However, we note a significant frac- 
tion of detections and lower limits higher than 300 keV: 
6 sources out of 20 have have Ec >300 keV. We conclude 
that a photoelectric cut-off in the 100-300 keV band is not 
an ubiquitous property of Seyfert galaxies. 



4.4. The iron line 

Equivalent width: We detect an iron Ka line in all the 
observations but one (NGC 7679). The equivalent width, 
EW, is in the range 100-300 eV for sources with low ab- 
sorbing column density (Nh < 10 23 cm -2 ), in agreement 
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Fig. 4. Observed equivalent width of the iron Ka line versus 
the absorbing column density. The two lines are the average 
equivalent widths for the subsamples with Nh < 3 x 10 23 cm -2 
and Nh > 3 x 10 23 cm' 2 , respectively. The shaded regions are 
delimited by the 1 a errors on the averages (for the low-iVn 
objects we excluded the point in the upper-left of the plot 
from the calculation. 



with previous studies, and with the values found in Seyfert 
Is (Gilli et al. 1999). In sources with a larger Nh, the con- 
tinuum at the energy of the line is partially absorbed, and 
the EW of the iron line is therefore expected to be higher. 
In Fig. H where we plot the EW of the iron line versus 
Nh- A correlation is present: sources with Nh > 3 x 10 23 
cm~ 2 have on average an EW higher that sources with 
Nh < 3 x 10 23 , even if the errors are too large to draw 
solid conclusions. 

This correlation suggests that the iron line is produced, 
at least in part, by a reflector different from the accretion 
disk. Otherwise, we would not expect any increase of EW 
with Nh, since the line would be as much absorbed as 
the continuum component. In Fig. ^ we plot the line EW 
versus the cold reflection ratio, R. A correlation is non 
clear, even if the errors are again too large to rule out an 
average increase of EW with R. We also note that in this 
kind of correlations, where the expected variation of EW 
is a factor of no more than 2-3, an important role could be 
played be the inclination of the accretion disk with respect 
of the line of sight: the primary continuum emission from 
an edge-on disk is weaker than that of a face-on disk with 
the same intrinsic luminosity. Therefore, if the iron line 
is partly produced by an outer reflector, the EW depends 
on the inclination angle. The effect is a further dilution of 
the correlations we are trying to investigate. 

Line energy and width: In all cases the measured 
peak energy is compatible with neutral or little ionized 
iron. In 9 cases out of 31, we measure a line width signif- 
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Fig. 5. Equivalent width of the iron Ka line versus the cold 
reflection component (normalized to the primary power law). 



icantly larger than 0. For all these cases we analyzed the 
model residuals to look for asymmetries in the line profile. 
We did not find any significant asymmetry, except than in 
one case, NGC 6300, which is discussed in the Appendix, 
and represents a convincing example of a relativistic line 
in Seyfert 2s. 

4.5. Soft excess 

In this work the soft emission is not analyzed in detail. A 
single thermal component, plus the warm reflection com- 
ponent, are used to fit the points at energies lower than 
the photoelectric cut-off. This gives on average good fits. 
The residuals at low energies are in no case high enough 
to alter the best fit values of the other high energy param- 
eters which we are more interested to. In the simple model 
A we fit the soft component with a power-law with a free 
photon index. The ratio between the absorbed and the 
unabsorbed power-laws, reported in Column 6 of Table 3, 
shows that the normalization at 1 keV of the "soft" com- 
ponent is typically a few percent of that of the absorbed 
component. This confirms the separation of our fits in two 
almost independent parts: the high energy components de- 
pends only on the data at energies higher than the cut-off, 
while the soft component is entirely determined by the low 
energy counts. 

5. Conclusions 

We presented a homogeneous analysis of 31 spectra of 20 
Compton thin Seyfert 2s observed with the BeppoSAX 
satellite. We fitted the spectra with multi-component 
models, including a cold and warm reflection, an iron line, 
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a low-energy thermal emission, and a primary continuum 
modeled with an absorbed power law with a high-energy 
cut-off. We also investigated the main correlations among 
the best fit parameters, in order to test the unified models. 
The main conclusions are the following. 

1. Testing Unified models: The spectra of the objects 
in our sample are in good agreement with the stan- 
dard unified model of AGNs: the primary continuum 
is well represented by a power law with photon index 
T = 1.79 ± 0.01, with a dispersion a(T) = 0.23. This is 
in agreement with previous studies on the hard X-ray 
emission of Seyfert 2s (Turner et al. 1997, Nandra & 
Pounds 1989, Smith & Done 1997, Turner & Pounds 
1989), and with the hard X-ray spectra of Seyfert Is. 
All the objects but one show significant cold absorp- 
tion by gas with column density Na > 10 22 cm -2 , in 
agreement with the optical classification of type 2 or 
1.8/1.9. The only exception, NGC 7679, is an intrigu- 
ing object, with no significant absorption in excess to 
Galactic in the X-rays, but also with no evidence of 
broad lines in the optical (see the Appendix and Delia 
Ceca et al. 2001 for further discussion). We also fitted 
the hard spectrum of our objects with a simple ab- 
sorbed power law, in order to compute the average ob- 
served photon index, r c ff. We found r o ff = 1.76 ±0.01, 
with a dispersion a(T e g) = 0.21. 

An iron Ka emission line at energies ~ 6.4 — 6.6 keV is 
found in all the objects but one (NGC 7679). In agree- 
ment with the unified model, the average equivalent 
width of the line is higher in objects with higher ab- 
sorbing column density. In 6 cases we found a broad 
profile of the line. In one case (NGC 6300) this pro- 
file is best fitted with a relativistic model. In the other 
5 cases, the statistics is not enough to test models of 
relativistic lines. 

2. A reflected component is present in most of the objects 
(17 out of 21). A debated issue on X-ray spectra of 
Seyferts is whether the reflection is due to the accretion 
disk or to a farther reflector. Comparing the variations 
of the direct and reflected components in the 7 sources 
with multiple observations, we find that the variability 
data are better reproduced assuming that the reflector 
is located far (> 1 light year) from the primary source. 

3. The ratio between the reflected and direct components 
is on average too high for a homogeneous reflector 
with the same column density as measured in absorp- 
tion in our spectra. This is an indication that the cir- 
cumnuclear medium is not homogeneous and proba- 
bly a Compton thick gas cover a significant fraction 
of the solid angle. This structure is also suggested by 
X-ray variability studies of Seyfert galaxies (Risaliti et 
al. 2002) and will be further tested in a forthcoming 
detailed study of variability of this sample of AGNs 
(Risaliti et al. 2002, in preparation). 
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Appendix A: Analysis of single sources 

NGC 526a: The same set of BeppoSAX data on this source 
have been analyzed by Landi et al. (2001). The results ob- 
tained by these authors are in agreement with ours, within the 
errors. In particular, it is confirmed the relatively flat contin- 
uum spectrum: the best fit values for the photon index of the 
intrinsic power law are in the 1.47-1.6 range, in agreement with 
our result, F = 1.54±g;g|. 

NGC 1365: NGC 1365 is one of the few sources for which 
our models do not provide a good fit (x?=l-26 in model C). A 
more careful analysis of this object was performed by Risaliti et 
al. (2000). According to this work, the PDS data of NGC 1365 
could be contaminated by the nearby active galaxy NGC 1386. 
Moreover, NGC 1365 shows an high and complex variability 
during the BeppoSAX observation. However, the best fit values 
in Risaliti et al. 2000 are in agreement with those obtained with 
our model C. 

NGC 2992: This source shows the strongest long-term vari- 
ability in our sample. Our analysis is in agreement with that 
of Gilli et al. (2000). 

IRAS 05189-2524: The BeppoSAX observation of this source 
was analyzed in detail by. Severgnini et al. 2000 The reflected 
component is not required from the fit (A% 2 = 2 with one more 
free parameter). The fit obtained with model B is in agree- 
ment with that of Severgnini et al. (2000) . The poor statistics 
(x? = l-2) is mainly due to an excess in the PDS data around 20 
keV. The statistics is however not good enough to understand 
whether this excess is real or not. We note that given the high 
noise in the PDS data, Severgnini et al. 2000 chose not to use 
these data in their analysis. 

NGC 4258: NGC 4258 is a nearby galaxy (estimated dis- 
tance ~ 7 Mpc) hosting a weak AGN. In this object the in- 
tegrated X-ray emission of galactic sources is not negligible 
with respect to the AGN contribution. Therefore, in models B 
and C we added an extra bremmshtralung component to take 
into account for the galactic contribution. The best fit temper- 
ature (kT=7 keV) and normalization are in agreement with 
those estimated in 5 different measures performed with ASCA 
(Reynolds et al. 2000). Our results is in agreement with the 
more detailed analysis of the same data performed by Fiore et 
al. (2001). 



NGC 6300: In Fig. Al we plot the residuals of the best fit of 



NGC 6300, obtained with model C, but without the Gaussian 
component that fits the iron Ka line. The asymmetry of the 
line is apparent: a clear "red wing" extends down to ~ 4 keV, 
while a broad blue wing is absent. We fitted this excess with 
a relativistic line model. We used the DISKLINE model in 
XSPEC, leaving all the parameters free. The best fit is obtained 
with a strongly relativistic profile: the best fit inner radius of 
the emitting region is 6.5 Shwarzchild radii; the inclination is 
~ 30°. The improvement of the fit, with respect to the simple 
gaussian fit, is A\ 2 = 27. We also tried a diskline modelem 
for a rotating disk, but the fit is worse, since this model can 
well reproduce the broad red wing, but not the narrow peak 
at E=6.1 keV. We note that a possible warning on this result 
could come from the high Nh of this source (Nh ~ 2 x 10 23 
cm -2 ). This implies that the continuum at the low energy peak 
of the line is partially absorbed. This could introduce system- 
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Fig. A.l. Residuals of the best fit model for NGC 6300, after 
the removal of the relativistic line from the model. 



intrinsic continuum (which is completely covered by gas with 
Nh ~ 1.4 x 10 23 cm -2 , in agreement with our results), is also 
absorbed by a second screen with Nh ~ 1.6 x 10 24 cm" 2 , and 
a covering factor of ~ 60%. 

NGC 7679: This source is peculiar within our sample, for it 
is the only one with negligible X-ray absorption. The optical 
classification (pure Seyfert 2) suggests that the BLR emission 
is obscured by dust by at least 3-4 magnitudes . This, assuming 
a galactic dust-to-gas ratio, would imply an absorbing column 
density of at least 6-7 10 21 cm" 2 . We can explain this dis- 
crepancy in two ways: a) the dust-to-gas ratio could be much 
higher than galactic, or the dust could be associated to a warm 
absorber; b) there could be a "hole" in the absorber of the nu- 
clear source, so that the X-ray emission from the accretion 
disk is unabsorbed, but the broad lines, emitted by the larger 
Broad line Region, are obscured. This scenario is not unplau- 
sible, given the recent evidence of a complex structure in the 
Seyfert 2 absorbers (Risaliti et al 2002) . Both the fits with and 
without the reflected component are acceptable (see Table 4 
and 5). In our analysis the best reduced \ 2 i s obtained with 
a steep F ~ 2.1 power law plus a strong reflected component. 
However, in a recent study of the same data, Delia Ceca et al. 
(2001) propose a best fit model in agreement with the one in 
Table 4, without reflection, and with a photon index F ~ 1.75. 



atic errors, making the reality of the red wing of the line less 
certain than it appears in our fit. 

Another possible scenario is the one proposed by Guainazzi 
(2001): a good fit of the MECS spectrum can be obtained ac- 
cepting very high values of the ratio R between the reflected 
and intrinsic component (R~ 4). In this way the flat profile of 
the cold reflection spectrum can reproduce the excess between 
4 and 6 keV which has been interpreted above as a relativistic 
wing of the iron line. Such high values of R are possible, as 
suggested by Guainazzi (2001), if the intrinsic flux is variable, 
since the observed reflection component is produced by radia- 
tion emitted at a different time with respect to the observed 
primary emission. 

ESO 103-G35: The BeppoSAX data of this source have been 
analyzed by Wilkes et al. (2000) . Even if the models used in our 
work are slightly different, the basic results are in agreement 
with these authors. The only noticeable discrepancy is about 
the high energy cut-off measurements. In the first observation 
we find Ec = 84I33 4 keV. This is the lowest measured cut-off 
energy in our sample, but it is significantly higher than that 
found by Wilkes et al. 2000 (E c = 29 ± 10 keV). 
NGC 7172: There are two published works on the analysis 
of the same data (Akylas et al. 2001 and Dadina et al. 2001). 
Akylas et al. 2001 propose a model with no reflection. Their 
best fit parameters are in agreement with our model B (Table 
4). Dadina et al. 2001 find that the inclusion of a cold reflection 
slightly increase the goodness of the fit. The best fit slope of 
the continuum is significantly higher in this model Our results 
are in agreement with Dadina et al. 2001. Our best fit model C 
(Table 5) is only slightly better than model B, and predicts a 
steeper intrinsic continuum (F ~ 1.9), than model B (r ~ 1.7, 
in agreement with Akylas et al. 2001). 

NGC 7582: Turner et al. (2000) present a detailed analysis 
of the same BeppoSAX data. The results are compatible, but 
these authors prefer a different interpretation for the high en- 
ergy emission in excess of the simple powerlaw, which we fit 
with a strong cold reflection component. Arguing that the re- 
quired reflection component is too high, they suggest that the 
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Fig. A. 2. Spectra and residuals for the whole sample. Residuals are calculated using model A (absorbed power law plus a 
second power law to fit the soft component). 
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